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Abstract 

c-| ! We formulate the Schwinger-Dyson equations in the ladder approximation for 

2D induced quantum gravity with fermions using covariant gauges of harmonic type. 
It is shown that these equations can be formulated consistently in a gauge of Landau 
type (for negative cosmological constant). A numerical analysis of the equations hints 
towards the possibility of chiral symmetry breaking, depending on the value of the 
coupling constant. 
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A very useful tool for the study of dynamical symmetry breaking and dynamical 
mass generation in quantum field theory are the Schwinger- Dyson equations. The stan- 
dard approach, which consists in working with this infinite set of integral equations has 
been developed in the pioneering works Jl|, , where the example of quantum electrody- 
namics was considered. By investigating the truncated version of the Schwinger-Dyson 
equations (the ladder approximation), the possibility of chiral symmetry breaking and 
dynamical fermion mass generation in QED could be demonstrated |l], (for a re- 
view, see the proceedings 0). The critical coupling constant in the Landau gauge has 
been found also. However, if dynamical fermion mass generation is taking place, then 
the non-perturbative Ward-Takahashi identity of QED is not satisfied. As a result, 
the dynamical fermion mass and the critical coupling constant in QED are very much 
gauge- dependent (for a recent discussion in an arbitrary covariant gauge, see ||, |j and 
references therein). 

If one is interested in further developing the Schwinger-Dyson formalism then 
more complicated models to study such equations, like those of quantum gravity, are to 
be considered. An example of an investigation of this kind for the case of 4D Einstein 
gravity coupled to fermions on a flat background has been presented in ref. ||, where 
by means of numerical estimations the possibility of chiral symmetry breaking has again 
been shown. 

In the present work we will study the Schwinger-Dyson equations in 2D quantum 
gravity |7j with fermions, on a flat background. The covariant gauge with two gauge 
parameters will be chosen and chiral symmetry breaking in 2D quantum gravity will be 
investigated numerically. 

The action of the theory under discussion is given by 

S = Sg + Sf, 

l r -o , — /„! 



S g -—Jd 2 x^{R-R + k 

Sf = JcPxy/^giV-fDfl, (1) 

where R is the two-dimensional curvature, \l/ the 2D spinor, A the cosmological constant, 



and D^ = d^ — (z/4)ur° a\, c is the 2D covariant derivative for spinors, where a, 



ab 



(i/2) [70, 7b], oj bc ^ is spin-connection, and e a will denote the vierbein. 

In the standard approach to string theory 0, one can start from a pure-matter 
theory (as given by Sf) in an external gravitational field, integrate then over spinors 
(this is easy to do in the conformal gauge) and get finally as a result 2D induced gravity, 
S g , where 7 is then specified. 

Here we are going to employ a more traditional approach, in which we will start 
from the theory (|]) and use the background field method on the flat background 

Hence, 7 in ([!]) is some given constant and we do not integrate over spinors. Expanding 
Sf on the flat background and working in the momentum representation, one easily 
finds that the interaction Lagrangian has the following form 

1 
Lmt = -f{p) [2W}»» - H»P») + (%- - 7(A))] ^ip)hT[k), (3) 

what corresponds to the fermion-graviton vertex 

hapu = T {^VXaVu-v - VXuVva - V^Vap) ■ ( 4 ) 

The gauge-fixing action will be choosen in the following form 

S 9 f = -^-Jd 2 xV=9^ (V^ p - (3V p h) (VX P - (5V p h) , (5) 

where a and (3 are the gauge parameters. (For a discussion of 2D induced gravity in a 
covariant gauge of the harmonic type, see also || |§). 

The quadratic part of the total action S = S g + S g f on a flat background is found 
to be 

S {2) = ~y J ' d 2 xh^H, upa hP°, (6) 

where 

H^upa = ^ V u P + 7j6 (Vpa^^u + T^VpVff) + &VpvVp* n 

1 A 

+ 76 yn^W V V p + rjppV u V a + r] va VpV p + l q V p\/ (J y a ) + -/ Wff , (7) 



being £1 = 2((3/ot — 1), £2 = 1 — ft 2 /a and £3 = — 1/at. The graviton propagator is given 
by the inverse operator H~^ and can be found using the algorithm of refs. [HJ, which 
yield 



4 1 / [ a + (3(l-2(3)f a 

A » vpa+ a-{2(3-l) 2 \{{3-l) 2 {k 2 -m 2 ) k 2 



G (k)--^' 1 (k)---L 1 l [a-hp^-^j a 



2a M 1 X f (2/?-l)[a + /?(!- 2/?)] a\ 



k 2 -m 2 * vpa a - (2(3 - 1) 2 I (/? - l)(fc 2 - m 2 ) fc 2 . 



M^/ 



+ 

Here 



4 1 / (1-2/3) 2 o_' 

A a-(2/3-l) 2 I (F-m 2 ) k 2 



-L^-Lpcr. (o) 



kpk u T _ 1 , , __ fc^ft,, 



-^pi/ 'Ipv 1 2 ' -'-'pupa ,-. \*->pp*->vo > *->po*-'vp) j ± pv , 2 ' \ / 

_ k^k u k p k a -}l(t P + r P + r P + t P ) m 2 _ aA 

The exact spinor propagator has the following form 

S- 1 (p) = A(p)p-B(p 2 ) } (10) 

where A and B are some unknown functions. Now we have at hand all the Feynman 
diagram elements: the exact spinor propagator, the free spinor porpagator, Sq 1 ^) = p, 
the free graviton propagator (§) and the vertex (£|). 



The effective potential for the composite fields |TTJ in the ladder approximation 
[0, |] can be written as 

Ves = -1 Sp (inS^S - S^S + l) + V 2 , (11) 

where V2 corresponds to the two-particle irreducible vacuum diagram, which follows 
from the vertex 

V2 = 4/(4^/(i^ Tr [r(p-q,q)S(q)T( q -p,p)G(p)]. (12) 

The Schwinger- Dyson equations (in the ladder approximation) correspond to the 
minimum of the potential (TlT| ) 

»2 

S-\p) - S \p) = -ij jJLr^fap- q)S(q)T pa (p, q - p)G^»°(p - q). (13) 



Using the explicit form of the spinor and graviton propagators flTU| ) , (|8[) , and the vertex 
(f|), one can get — after performing Wick's rotation and the angular integration (we 
drop the details of these straightforward but very tedious calculations) 



+9 



N f M 2 


{lo dX 


\\n(i 2 (x) 1 B2{X) 


:ft — „ 

8n 


\ X 


dx 


f 1 d y 



v; 



o xA 2 (x) + B 2 (x) Jo yA 2 (y) + B 2 (y) 



A(x)(A(x) - l)x + B 2 (x) 
xA 2 (x) + B 2 {x) 

[A{x)A{y)K A {x, y) + B(x)B(y)K B {x 



(14) 



where Nf is the dimension of the fermion representation, M the momentum cutoff, 
x = p 2 /M 2 , y = q 2 /M 2 , A(x) = A(p 2 ), B(x) = B(p 2 )/M, and g = 7 /(647r) and 
I = A/(2M 2 ). The explicit expressions for Ka and Kb are very complicated for arbitrary 
a and /3. Moreover, if a ^ the Schwinger- Dyson equations contain the infrared 
divergences caused by the graviton zero momentum. Let us give some examples of Ka 
and Kb for different choices of the gauge parameters. 



1. Gauge with a arbitrary, /3 = 1/2. 



K A (x,y) = -{(l-3a)(x + y) + 



a[7{x 2 + y 2 ) + lOxy + 3al(x + y)\ + A(x + y)(x- yf jl 



K B (x,y) = 5a - 1 



'x + y + al) 2 — Axy 

x 2 + y 2 + Qxy + 4(x + y)(x — y) 2 /l 1 

i i i ' 

\x-y\ J 

5a[2(x + y) + al] + 4(x - y) 2 /l 2(x + y) + 4(x - y) 2 /l 



'-. (15) 



(x + y + al) 2 — Axy \ x V\ 

Here A > and I = A/(2M 2 ), and in this gauge one finds infrared divergences in the 
Schwinger-Dyson equations (i.e., at the lower limit of the integrals in (|I4|)). 

2. Let us now consider a gauge of Landau type (a = 0, (3 arbitrary), in which the 
Schwinger-Dyson equations do not contain infrared divergences. There one finds 



K A {x,y) 



ff 2 



+ 



2 K/3-1) 2 

4(4/3 - 1; 
W - i) 



[x + y) 



(5[A(3xy + 2(2/3 - l)(x - y) 2 + /3(x + y)(x + y + ^ 2 )] 



(x + y)(a;-y) 2 +f 4 ( ^_ 1 i ) ) j (ar-y)* 



(/^-l) 5 



[x 2 + y + /i ) - 4xy 



-1/2 



8/3 



K B (x,y) 



1 



09 -1) 



Z(2/? - 1 

-[2(x + y) + fi 2 } + 



(x + y)\x-y\ + 



>^ , 4(x-y) 2 



/(2/3-1; 

4|x — y\ 



4(/3 - 1) 
1(2(3-1) 

x 2 + y 2 + ^ 2 ) 2 
P 2 



x-y\ 



4xy 



-1/2 



1(2/3-1) (P-1) 2 ' (16) 

Here A < 0, / = -A/(2M 2 ) and /i 2 = (4/1) [(2(3 - l)/((3 - l)] 2 . Notice that for a = 0, 
/3 = 1/2, the theory contains again infrared divergences, because in this case 



K A (x,y) 
K B (x,y) 



1 ( x 2 + y 2 + 6xy" 
- [x + y : : 

2 V \x-y\ 



x + y 

\x-y\ 



1. 



(17) 



Observe also that, in principle, one expect to find more complicated covariant gauges 
which are free of infrared problems in the region where the cosmological constant is 
positive. 

Starting from eqs. ( |T3| ) and integrating over the angles one can show that the 
functions A and B must obey integral equations of the following form 

A(y) 1 



A(x) 
B(x) 



1 + g dy 
Jo 



yA 2 (y) + B 2 (y)x 

B(y) 



K A (x,y), 



9 lo dy yA 2 (y)7B*(y) KB ^ y) - 



It is not possible to solve these equations analytically (We will discuss here the case 



of the physical Landau- type gauge (16) only, where no IR divergences appear in the 
theory, in order to avoid the introduction of any IR cutoff). We present the result of a 
numerical calculation, obtained by using an iterative procedure (in close analogy with 
H). We consider two types of trial functions 

(a) A°(x) = d, B°(x) = 0, 

(b) A°(x) = ci, B°(x) = c 2 , 

where c\ and c 2 are some constants between and 1. We will also fix the values of g, 
I and (3. The functions A°(x) and B°(x) can then be taken as the starting point of a 



self-consistent iterative calculation of the form 

A* +1 (x) = 1 + g dy [ ^ , -K A (x,y), 

Jo yA l (y) + B* (y) x 

B i+1 (x) = gfdy .. 2 *'i y L, K B (x,y). (19) 

Jo yA l (y) + B l (y) 

The sequences formed by the {A l (x)} and {B l (x)} are expected to converge towards 
the functions A(x) and B(x), respectively, which are the sought for solutions of fll8f ). In 
practice one can judge the degree of convergence of these series by the smallness of the 
squared norms of the differences A l+1 — A 1 and B' l+1 — B l , which we set at 10 -4 — 10 -6 
in our calculation. If, for the given g and /, there are solutions of both types, (a) and 
(b), only the most stable of both by V c s fllPf) is to be chosen as the one corresponding 
to the true vacuum. 

We have executed this algorithm to solve (|T9"D, starting from the trial functions 
(a) and (b), for fixed I = 4, j3 = 1/3 and varying g. For very small g's, both types 
lead to curves close to A(x) = 1, B(x) = 0, i.e. the chiral symmetric solution, as was 
to be expected. As g increases, the value of V e ff for the chiral solution of symmetric 
type (a) appears to be slightly higher than the corresponding one for the non-symmetric 
solution. In particular, for g = 0.1 (see Fig. 1), the chiral symmetric solution is the 
preferred one. For g = 0.2 or g = 0.3 (see Fig. 1 again, where typical curves for A and 
B are presented), one can see that the chiral non-symmetric solutions are preferable. 
Hence, we see clearly that the Schwinger- Dyson equations for 2D gravity with fermions 
may have chiral symmetry breaking regimes in the covariant gauges. 

We will now say a few words about the regime of the Schwinger-Dyson equa- 
tions corresponding to a theory with positive cosmological constant. In this case, when 
working in the covariant gauges under discussion, one encounters problems related with 
infrared divergences. We will use the conformal gauge 

g^u = e v r]^ u . (20) 

In this gauge one finds the Schwinger-Dyson equations (|T3"D and ( |l"8D with the functions 



Ka and Kb |T2| gi yen by 



Axy + (x + y)(x + y + l-J(x + y + l) 2 -Axy) 

K A (x,y) = , 

2y (x + y + l) 2 - Axy 

2(x + y)+l- J(x + y + l) 2 - Axy 

Kb(x,v) = Vl -■ (21) 

(x + y + I) 2 — Axy 



Numerical solutions of the corresponding Schwinger- Dyson equations can be obtained 



as above (see ref. [0). Typical curves for I = 0.5 and varying g are shown in Fig. 2. 
Here we observe again the possibility of chiral symmetry breaking. 

Summing up, we have studied the Schwinger-Dyson equations corresponding to 
2D gravity coupled with fermions in a covariant (harmonic) gauge. Numerical anal- 
ysis of the equations show clearly the possibility of chiral symmetry breaking in the 
region with negative cosmo logical constant, where the Schwinger-Dyson equations can 
be consistently formulated in a gauge of Landau type and no infrared divergences ap- 
pear. In the region of positive cosmological constant, the analysis done in the conformal 
gauge shows as well the possibility of chiral symmetry breaking. The results of the 
numerical analysis of the solutions (and the Schwinger-Dyson equations themselves) are 
certainly gauge dependent. Currently there is no way to solve such a drawback of the 
Schwinger-Dyson equations, even in the case of renormalizable theories, as QED, where 
the Ward-Takahashi identities have a quite simple form. We have nothing to add here 
that can help to resolve this general problem of gauge dependence [|, |]. Our purpose 
has been simply to show that chiral symmetry breaking is indeed possible in 2D grav- 
ity theories with fermions in different gauges. That this is actually the case has been 
realized by means of a rather straightforward numerical analysis of the corresponding 
Schwinger-Dyson equations. 
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Figure captions 

Fig. 1. Plot of the functions A and B obtained as the (a)-type solutions for g = 0.1 
and (b)-type solutions for g = 0.2 and g = 0.3 keeping / = 4 fixed. 

Fig. 2. Plot of the functions A and B obtained as the (a)-type and (b)-type solutions 
for g =0.1, 0.2 and 0.25, keeping I = 0.5 fixed. Notice how B deviates more and more 
from the g = solution (B(x) = 0) as g increases. Although not shown in the figure, 
the curve keeps going up for larger values of g. 
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